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Abstract 

We consider the evolution of a network of strings in an expanding universe, allowing for the 
formation of junctions between strings of different tensions. By explicitly including, in the velocity- 
dependent evolution equations for the network, kinematic constraints associated with the formation 
of Y-shaped string junctions, we show how they lead to scaling solutions in regimes where they 
would not otherwise be found, thereby extending the range of parameters which lead to scaling. 
By incorporating these constraints we are able to study their general behaviour for networks with 
cosmic superstring interaction rules, and predict the scaling densities expected by these networks. 
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I. INTRODUCTION 



Understanding how a network of cosmic strings evolves in our universe has been an 
area of interest for a number of decades following the pioneering work of Kibble 
(for a review see 0, J]). Being a complicated non- linear system, various approaches have 
been adopted including numerical simulations Q, B as well as analytic approximations Q, 
0, y, 9|. The interest in the results of these simulations wained when it became evident 
that cosmic strings by themselves were never going to provide the seed primordial density 
fluctuations [10], but recently interest in them has resurfaced thanks to the realisation 
that in a reasonabl y la rge class of Superstring motivated models it is possible to form a 



network of strings UJ, Q, [l4|, that can live as long as the age of the Universe, be 



of cosmological length and survive a period of early universe inflation. Moreover these 
cosmic superstrings could possibly be detected, thereby providing the first direct evidence 



of string theory in nature (for recent reviews of cosmic superstrings see 



m 



16 



12, 



y 



19|). This has spawned a renewed interest in the subject, in particular there is now a 



need to understand how a network of cosmic superstrings evolves. The fact that there 
can now be more than one type of string has added a new ingredient, the possibility of 
forming junctions when two different types intersect, the junction marking the point where 
a third string emerges. A number of approaches have already been adopted to study this 



problem. On the numerical side, network evolution with junctions has been studied in 20|, 



21 



22 



23j. Preliminary work based on non-abelian field theory models had suggested that 



networks with junctions rnay be obstructed from sca ling j2l| (network frustration), but, more 



recently, analytic 2J, |25|, |26| and numerical 22|, |23|, |27|] evidence has accumulated favouring 
a picture in which these networks^ reach scaling, while the various network components are 
allowed to interact with each other forming 3-string Y-type junctions. The development 
of these field-theory analogues of cosmic superstrings have involved either combinations of 
interacting abelian models 



23 



28 



291] or through non-abelian models 



20 



21 



22 



27|. In 



both cases the networks that form admit trilinear vertices, hence junctions where the usual 
intercommutation p.opertie^ of the strings no longer apply. ^ 

On the analytic front, the first successful, comprehensive model was that of [25|], which 



^ Similar results are also found in field theory simulations of domain wall networks with junctions 



30 



31 



3- 



2 



was based on the velocity-dependent string evolution models of Martins and Shellard 



In 



],Q. 



25l | the authors assign a single correlation length and velocity to all network components 



and follow the evolution of the number density of each of the different species, finding generic 
scaling under the assumption that the energy associated to the formation of junctions is 
rapidly radiated away and decouples from the string network. In order to make the basic 
model more realistic, in reference 26j, the authors extended the velocity-dependent models 



by associating a different correlation length and velocity to each network, and enforcing 
energy conservation at junction formation. While once again, scaling solutions are readily 
found in this model, and are in qualitative agreement with those of reference 25(], the 
scaling behaviour is seen not to be generic in the model of 2^. In fact, for scaling to occur, 
it requires that the terms describing the formation of junctions be subdominant - or at most 
comparable in magnitude - to the corresponding self-interaction terms. 

As well as understanding the dynamics of the string network, properly incorporating the 
kinematics of strings that can form junctions is equally important and this area has recently 



started receiving attention 



33 



34 



35l | . Studies of the collisions of Nambu-Goto strings with 



junctions at which three strings meet have shown that the exchange to form junctions cannot 
occur if the strings meet with very large relative velocity. For the case of non-abelian strings, 



rather than passing through one another they become stuck in an X configuration [3j, |36| , 
in each case the constraint depending on the angle at which the strings meet, on their rela- 
tive velocity, and on the ratios of the string tensions. Calculations of the average speed at 
which a junction moves along each of the three strings from which it is formed yield results 
consistent with the observation that junction dynamics may be such as to preferentially 
remove the heavy strings from the network leaving a network of predominantly light strings. 



In 



37( 1 the authors extended the analysis to include the formation of three-string junctions 



between {p, g)-cosmic superstrings, which required modifications of the Nambu-Goto equa- 
tions to take account of the additional requirements of flux conservation. Investigating the 
collisions between such strings they showed that kinematic constraints analogous to those 
found previously for collisions of Nambu-Goto strings apply here too. 

A number of these constraints that have emerged from analysing the modified Nambu- 



Goto equations have been checked through 
which can also lead to junction formation 



n numerica. 



38 



39 



40 



simulations of field-theory strings 



4l| . In most cases good agreement 



has been obtained with the analytical predictions, although the results of [23j show some 
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surprising features, and of course some important differences emerge due to the fact that 
field-theory strings have interactions between them that are not found in the Nambu-Goto 
case, hence some of the detailed numbers differ. 

In this paper we show how to consistently include the kinematic constraints associated 



with Y-shaped string junctions 



33 



34 



37l | into the model of [26|. A first step towards 
this goal has been taken in reference |42| by introducing a velocity cut-off beyond which 
colliding strings fail to exchange partners. Here, we take into account the dependence of the 
constraints on both the velocity and orientation of the colliding strings, and demonstrate how 
these constraints provide the required suppression factors for the relevant junction-formation 
terms, thereby facilitating scaling in regimes where scaling behaviour would not otherwise 
be observed. By incorporating these constraints in the velocity-dependent evolution models 
we are able to study their general behaviour for networks with cosmic superstring interaction 
rules. 

The rest of the paper is as follows. In section [III we introduce the velocity-dependent 
string evolution models, and in section IIIII we describe the kinematic constraints that apply 
to {p, g)-junction formation in cosmic superstrings. This is followed in section IIVI with the 
incorporation of the constraints in the string evolution models. Our key results are presented 
in section |V] before we conclude in section |Vll 



II. VELOCITY-DEPENDENT STRING EVOLUTION MODELS 



Let us begin by briefly reviewing the model of reference [26|, which will be the starting 
point of our present discussion. This is a phenomenological model for the quantitative 
description of the cosmological evolution of cosmic string networks that are composed of 
different string types, generally of unequal tensions, allowed to interact among each other 
forming trilinear vertices (Y-type junctions). Each network component i is assumed to 
be Brownian (that is of a random walk structure) and is characterised by its tension /ij, 
correlation length Lj and root-mean-squared (RMS) velocity Vi, the latter two quantities 
being functions of cosmological time t. Due to the Brownian structure of the networks, 
the correlation length Li{t) defines, at any time t, a network energy density component 
p,(t)=/i,/L2(t). 

Starting from the Nambu-Goto action in a Friedmann-Lemaitre-Robertson- Walker 



4 



(FLRW) background with scale factor a{t) and averaging over the string worldsheet, one 
can derive evolution equations for pi and Vi that account for string stretching and dilution 
as well as velocity redshiftin g d ue to cosmic expansion. One then introduces additional 
phenomenological terms , |26j describing the interaction among strings of the same type, 
producing loops, as well as cross interactions among different types of strings leading to the 
formation of Y-type junctions with new segments (links/zippers) connecting the colliding 



strings. The resulting equations are 



26|: 
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(2) 



where c, parametrises the self-interaction efficiency of strings of type i producing loops^ and 
dfl- = d^-^ denotes the efficiency parameter for the process in which strings of type i and j 
interact to produce a segment of type k. The mean velocity Vij is the average relative RMS 
velocity between strings of type i and j, and C-^jit) is the average length of links/zippers of 
type /c, produced by interactions between strings of types i and j around time t. In the 
second equation, ki is the curvature parameter, which indirectly encodes information about 
the small-scale structure on strings and can be expressed as a function of the velocity Viit). 



In the relativistic limit one finds 



9|: 



ki 



2^21 -8t 



(3) 



TT l + 8t;f 

Under the Brownian network assumption, the radii of curvature Ri are equal to the cor- 
responding correlation lengths Lj, so one should set Ri = Li in equation ([2]). Finally, the 
parameters b^j are introduced here in order to be able to switch between the model of ref- 

in which the energy liberated by the formation of junctions 



erence 



26| (where b' 



4)' 



25| 



is redistributed in the network as kinetic energy, and a model analogous to reference 
(corresponding to bfj = 0), where all this energy is assumed to be radiated away. 

One can use this general class of models in different contexts (Zjy-strings, non-abelian 
strings, cosmic superstrings) to describe networks with multiple string components that 



This can be expressed as an integral of an appropriate loop production function over all relevant loop 
sizes Isi. 
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interact with each other according to rules encoded in the coefficients dfj,Ci. The string 
tensions /i,, interaction rules (ifj,Cj and link/zipper length are model dependent. 

cosmic superstrings, for example, the interaction rules are well- 



In the case o 
understood 



25 



43 



44| . There are three possible outcomes of the crossing of a {p, q) and a 



{p'j q') string: they can either pass through one other or zip in two different ways, producing 
either a {p+p', q + q') or a {p — p', q — q') zipper segment. Given that the strings have zipped, 
the probability that the additive/subtractive channel is followed is given by: 

2 \ ^\ {p^gi + g2)i/2 + ^/2)i/2 y y ' y ) 

where we have assumed that the Ramond-Ramond scalar present in the analysis of 43|] is 
zero. Therefore, the coefficients d^^ in equation ([T]) can be expressed for the case of cosmic 
superstrings as: 

l(j>±p' ,q±q') _ 7 p± (r\ 

"'(p,q),(p',q') - "{p,g),{p',9')Mp,g),(p',g')' y^ 

where the factor (i(p^q)^(p'^g/) corresponds to the probability that the crossing is non-trivial 
(i.e. the strings do not pass through one other but rather form a zipper following either 
of the two above mentioned channels). This factor depends on the compactification and 
string coupling Qs, and for different string types lies in different ranges, which have been 



investigated in detail in references 



43, 



441 ]. We will treat these factors as free parameters 



chosen in the appropriate ranges, which typically lie between 10" < c?(p,g),(p',ij') < 1. 

The string tensions are also known. For example, for a string carrying charges {pi , g^), in 
flat spacetime, its tension is given by: 



n^ = —^phl + <i (6) 

9s 

where fip is the tension of the fundamental (1,0) string (F-string) and, again, we have 
assumed the Ramond-Ramond scalar to be zero. 

To close the system one also needs to specify a functional form for the average length of 
the produced links/zippers d^it). For superstrings the formation of junctions happens by 
zipping along the length of colliding strings, so the natural choice is to take (-^{t) equal to 
the smallest of the two correlation lengths Li{t), Lj(t). As in [26|, we will use 

LjLj 
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a simple expression returning a value smaller than (but close to) the smallest of the two 
correlation lengths. Other choices can be made, depending on the model, but as long as 
£ij{t) scales as some mean of the correlation lengths, these choices do not significantly alter 
the results 26|], especially in the case when the junction formation terms are subdominant, 



as will be the case here. 

Equations (JIHT]), truncated at any finite number of string species, can be solved numer- 
ically for any choice of the free parameters Cj, d^j. We now proceed to describe the micro- 
physical kinematic constraints for Y-junction formation, which we will incorporate into our 
string evolution models in section IIVI As we will see, including these constraints in our 
models will lead to a systematic suppression of the junction formation parameters d'^j. 



III. KINEMATIC CONSTRAINTS FOR {p, g)-JUNCTION FORMATION 



In [37j the authors extended arguments originally presented in j33| and [3J] to consider 
the kinematic constraints that apply to the formation of [p, g)-junctions. As these results 
are vital for our analysis we briefly review them here. Considering the collision of straight 
strings with charges (^1,^1) and {p2,(l2) meeting at an angle a and traveling with equal 
and opposite velocities v, when the strings collide, they may become linked by a string with 
charges (ps, qs) = — (pi +P2, Q'l + Q'2)- Imposing the requirement that the length of the joining 



string must increase in time leads to the strong kinematic constraints 



37j: 



where 



f^{v, a) =A,{1 + v'Y + ^2(1 + v') + A3 < 



(8) 



^3 
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7^2 -2 



-2 



COS a — sm a — /i_ cos a 



0-2 -2 2 -4 /o 2 i\-2 -2 

2jj,^j^_ COS a — jj,^ — [2 COS a — lj/i_,_/i3 



— 4 —2 —2 
/i3 - 



(9) 



and fi± = fii± Ji2- 

The inequality ([8]) leads to a condition on the values of v for which junction formation is 
possible, namely: 



< < vl{a), 



(10) 



7 



where the critical velocity, Vc, also depends on Jii and /i2, hence on the charges of the two 
colliding strings. In j34j], it was shown that ti™^^ < 1, (where v^^^ is defined as the maximum 
value taken by Vc as a varies from to 7r/2), only if 

fil>fi+\fi-\ = \fil-fil\. (11) 

Hence, if this condition is satisfied by the tensions of the joining strings, then there is a 
velocity v^^^ above which the two colliding strings will pass through each other rather than 
forming a junction. 



A couple of examples may help illuminate the effect of the constraints [37|]. For the 
collision of equal tension strings (for example (1,1) or (1,-1) strings), then /i_ = and the 
constraint Eqn. [8] becomes 



Vl-v^cosa > (12) 

Equation (fT2l) holds as long as the triangle inequalities are satisfied amongst the three 
tensions, namely that for all i 

Vi > (13) 

where, for example, 

1^1 = ^2 + fis - fil (14) 

(and cyclic permutations). For gs ^ 1, the largest region of the {a,v) plane is available to 
the collision of (^1,^1) with {pi,—qi) strings which produce a joining string with charges 
(— 2]9i,0). Such a string is a p-string which is light, indicating that it is much more likely 
to form a light joining p-string than a heavy joining g-string. There is a symmetry present 
which means that if the second string is of the {—pi, qi) type, then because of the invariance 
under p —>■ q and Qg —>■ I/qs (basically because F and D-strings are symmetric under string 
theory S-duality where weak coupling is replaced by strong coupling, i.e. gs ^/ds), then 
the linking string is a D-string. 

A second example is the collision of an F, (l,0)-string with a D, (0,l)-string. Such an 
event is important because it forms the building blocks for general (p,q) string collisions. Of 
course the third string formed at the junction is a (1,1) string. 

In js^], under the assumption that an x-link forms and 2 cos^ a > 1, the condition Eqn. ([S]) 
for junction formation is shown to become < f ^ < where 



2 (1 + 5ff) - 4cos^asin^a(l + 5fs)^ + y(l +5f2)2 _ 4cos2asin^a(l - 5f2)2 
^ 2cos2a(l + c/,)2(2cos2a- 1) ' 
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One of the interesting features that emerge is f c = when a = 7r/4, independent of Qs- For 
a > 7r/4 no a;-hnk can be formed in this case, whatever the values of Qs or v. 

As Qs 0, in which case f c = 1 then half the [a, v) plane is allowed implying physically, 
that a very heavy D-string can always exchange partners with a light F-string, with the 
third string moving with a velocity approximately equal and in a direction almost parallel 
to the incoming heavy D-string. A second limit is when Qs ^ 1 in which case the F-string is 
almost as heavy as the D-string and /i_ ^ 0, which corresponds to the case described above 
for equal tension strings. 

It is clear from this brief discussion that the presence of the kinematic constraints has a 
big effect on the ability of colliding strings to form junctions. We will now incorporate these 
constraints in the string evolution models and see how they have a dramatic effect on the 
scaling solutions obtained. 



IV. INCORPORATING THE CONSTRAINTS IN OUR STRING EVOLUTION 
MODELS 



The constraints described in the previous section act as to "switch off" certain interaction 
processes with parameters {v, a) outside the kinematically allowed range. These constraints 
have been derived from the Nambu-Goto equations and so they only (strictly) hold in the 
thin-string approximation. However, they have recently been found to be in excellent agree- 
ment with field theory simulations, based on Abelian-Higgs models 39|, |40|, although of 
course there are specific differences arising from the particular nature of the strings and 
their interactions (see also [23]). One can think of these constraints as window functions 
which, for each interaction process involving certain types of strings, restrict the allowed 
region of {v, Q;)-parameter space in which junctions can form. Thus, in a phenomenological 
approach like that of section [Tll where only the average effect of string collisions over the 
whole parameter range is considered, one should be able to integrate these constraints over 



the full range of the parameters {v, a) SJ, |45| . This would effectively introduce a suppression 



in the relevant coefficients of equations ([I])- ([2]) that would (in general) be different for each 
zipping process, since the constraint curves depend on the tensions of the strings involved. 
The aim of the present section is precisely to incorporate these constraints in the network 
evolution models of section [ITl In particular we will integrate the constraints of section IIIIl 
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modelled as window functions, over an appropriately normalised probability distribution in 
(f , Q;)-space. The result will be a number less than unity, since what we are doing is simply 
eliminating that part of parameter space which does not satisfy the constraint. In other 
words, we are effectively switching off a number of kinematically forbidden interactions that 
would otherwise be included in the models. As the string network is to a good approximation 
Brownian, the strings are randomly oriented so we will assume a fiat probability distribution 
for the collision angles a. On the other hand, the velocities v - for scaling networks - are 
peaked at a particular value f^, with a Gaussian probability distribution [7|, |46|. The vari- 
ance (Tt, depends on the tensions of the 3 strings involved (in fact on /ii and fl = ^2 — fJ'sY, 
but lies in a narrow range, for example crj G [0.239, 0.275] for 0.2 < fii < 1.4, p, < /ii. Thus, 
in the following we will take a universal variance of crj = 0.25 for all interactions. This can 
be easily generalised to o"„(/2), but the results will not depend strongly on this choice. 

Therefore, if ffi{v, a) is the constraint function in equation ([HD, so that the kinematically 
allowed region is given by /^(f , a) < 0, then the relevant suppression factor for the process 
where strings of types i and j interact to form a type k segment will be 

S^, = \ \l &{-Mv, a)) exp[{v - v,,y/a^]dadv < 1 , (16) 

where 9(x) is the Heaviside (Unit Step) function and Vij is the RMS relative velocity between 
strings of types i and j, that is 



= v/([v.(t)-v,(t)]2) = ^v,{ty + v,{ty . (17) 

Thus, the suppression factor (|T6|) depends on the particular interaction process through 
the tension dependence of ffi{v,a) and the relative velocity Vij"^. As shown in Fig. [1] these 
suppression factors are numerically very different for different interactions, so it is important 
to take them into account in macroscopic evolution models. 

The constraints can be then implemented in the (multi-)network evolution models of 
section [IT] by the replacement: 

(18) 



^ See equation (64) of |34| for (x^) = f{ij.i,fj.). One can do the same calculation for (Vii^)^, though in this 

case the result cannot be expressed analytically. 
^ Also through the variance if the tension dependence of is taken into account. 
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FIG. 1: The allowed parameter range in (w, a;)-space for different interactions, obtained from 
Eq. ([8]). The interactions shown are (starting from top left): (0,1) + (1,1) (1,0), (1,0) + (1,1) 
(0,1), (1,0) + (2,1) ^ (1,1), (0,1) + (1,1) ^ (1,2), (1,0) + (1,1) ^ (2,1) and(l,0) + (2,l) ^ (3,1). 
Upper plots correspond to the subtractive channel (making lighter string states) while lower ones 
to the additive (making heavier strings). The area of the allowed region depends strongly on the 
particular interaction, the subtractive channel (upper plots) being favoured over the additive one 
(lower plots). 

Clearly, with this replacement, the coefficients of the 3-string interaction terms in equations 
(II])-(I2D depend non-trivially on the unknown functions Vi{t) through the "parameter" Vij 
(see equations f|T6|l and f|T7j) ). so the system is numerically heavy to integrate. However, for 
scaling networks, one can solve the system iteratively, treating Vij S cLS constant parameters 
in each iteration. In the ffist step, one solves the system with initial guess-values assigned 
to Vij^s, and then solves the system again, replacing the initial guesses for Vij's with the 
corresponding values obtained by inserting the returned scaling values of fi(t)'s into equation 
(IT7|) . Repeating this procedure (in practice two or three iterations are needed), one can 
match the input Vij's to those obtained at scaling, to the required accuracy. In the next 
section we follow the above procedure to solve the system in the case of cosmic superstrings 
for different sets of values of Cj's and di/s. 
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V. NETWORK RESULTS 



The general behaviour of multi-tension string networks - without taking into account 



the above kinematic constraints 
contexts 



20 



21 



22 



231, 



25 



26 



27 



las been investigated by a number of authors in different 
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2^ and in particular, networks modelled by equations 



([TJ-Q have been studied in [26[. Here we will concentrate on the effects of the kinematic 
constraints of section IIIH which affect the network dynamics through the fact that they 
rescale the d^j terms (ITSl) . 

Let us first consider a simple model in which the interaction terms in the velocity evolution 
equations ([T]) are switched off, bfj = 0. Physically, this corresponds to the assumption that, 
during the binding process between strings of type i and j, forming a (lighter) string segment 
of type k, the energy difference arising from the disparity in string tensions is efficiently 
radiated away. This may be plausible for cosmic superstring networks, where binding is 
accompanied by the production of microscopic F-strings that can damp energy away. It was 
pointed out in reference |25i] that this mechanism alone can be sufficient to drive the network 
to a scaling regime, that is, with the assumption that string binding dissipates all the energy 
imbalance associated with junction formation processes, such junction-forming multi-tension 
networks can reach a scaling regime even when self-interactions among segments of the same 
type (which for single string networks provide the dominant energy-loss mechanism through 
loop production [3l]) are ignored. 

The system ([I])-([2]) with the interaction rules (|1D - ([5]), and b^j = generalises the model 



of 



25l | in that it assigns a different correlation length and a different RMS velocity for each 



of the network components, as described in section [TTl Even in this more general class of 



models the main result of 



25l ] remains true and scaling is generically reached, with heavier 



string components being systematically less populated. This is because the factors Pl^ in 
the coefficients d^j of ([T]) favour the subtractive over the additive interaction channel, so 
lighter strings acquire larger number densities compared to heavier ones. However, the 
relative abundance among strings of different types depends quantitatively on the numerical 
values of Cj - parametrising loop production - and dij = d\^/ - corresponding to binding 
interactions. Therefore, the kinematic constraints on string junctions, described in the 
previous section, can significantly affect the relative abundances of the most populated, 
lighter network components through the suppression factors f ll6l) . f[T8l) . which should therefore 
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be included in analytic models of this type. 

Figure [2] shows the scaling densities of such a network with string coupling = 0.3, 
for Ci = 0.23^, dij = 1^ (top panel) and for Cj = 0.023, dij = 0.1 (bottom panel), in the 
radiation era. Although the results here are for the radiation era, similar conclusions hold 
in the matter era but with different specific values for the scaling densities and parameters 
Ci and dij. The formally infinite set of equations ([I])-([2]) has been truncated at 4-string 
composites, leaving only seven species of distinct tension, namely (1,0) or F-strings, (0,1) 
or D-strings, (±1,1), (2, ±1), (±1,2), (3, ±1) and (±1,3), shown in black, blue, red, green, 
orange, pink and yellow respectively. The plots on the left have been obtained without 
taking into account the junction constraints of sections lllll and HV] while those on the right 
demonstrate the effect including these constraints have on the string scaling densities. Note 
on the left panel that, for these models (6^^- = 0), scaling is reached even though the cross- 
interaction coefficients dij are significantly larger than the self-interaction ones. Also note, 
that if we turn off the parameter determining the efficiency of cross-interactions, (i.e. we 
set dij = 0), then a network forms which is a tangle of different string species with tensions 
given by ([6]) that do not interact with one another. Such a network would reach a scaling 
regime in which all network components independently scale due to their self-interactions 
parametrized by Cj. For q = 0.23 (all i) all network components would asymptotically reach 
a normalised string density of = t^pi/^i ~ 14 in the radiation era, while for Cj = 0.023 
the corresponding value would be A^ ~ 980''. 

The overall reduction of scaling densities in the left panel with respect to the correspond- 
ing values for dij = just mentioned, reflects the impact of choosing values of dij > Ci. 
In that case the dynamics is dominated by energy losses associated with junction forming 
binding processes. If we also assume that the energy can be radiated away, then the energy 
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This value corresponds to that obtained by fitting Nambu-Goto simulations of single string networks 



^ The dij^s are chosen to be significantly larger than the Ci's to demonstrate that, for our b^j — model 
which is analogous to that of reference [2q . scaling can be found even though the junction formation 
terms are dominant. This is to be contrasted to the energy-conserving model b^j = d*^-, where scaling 
cannot be achieved for such large values of dij . 

^ This enhancement of string density resulting from a suppression of string self-interactions has generated 
significant interest in networks of cosmic superstrings, because their intercommutation probabilities can 



be suppressed compared to field theory strings ll|, |43| | . This in effect reduces Ci so cosmic superstring net 
works may . 
simulations 



works may have much larger number densities than their field theoretic counterparts. However, numerical 
suggest that this effect is not as strong as originally anticipated. 
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FIG. 2: Evolution of normalised string density N = in the radiation era for the lightest 

components of a network described by equations (HI),©,®,© with h^^ = and for Qs = 0.3. 
F-strings are shown in black, D-strings in blue, (l,ibl) strings in red, (2, ±1) in green, (l,ib2) 
in orange, (3, ±1) in pink and (l,ib3) in yellow. Plots on the left do not take into account the 
effects of junction constraints, while in the plots on the right, these effects are included through the 
coefficient rescalings p^ . Upper plots are for q = 0.23, dij = 1, while lower plots for q = 0.023, 
dij = 0.1. 



loss due to bindin g ca n be greater than the corresponding losses to loop-production. As 



noted in reference 



26| - and as will be demonstrated below - this is not observed in the 



case of energy conserving models, h^j = dfj, where the energy liberated by forming junctions 
is redistributed in the network through momentum exchange: such networks simply fail to 
reach scaling for large dij. The important effect of cross-interactions dij ^ 0, subject to 
quantitative modifications due to the inclusion of the junction constraints (right panel), is 
to suppress the abundance of heavy species relative to light ones through the 'selection rule' 
(jl]), which strongly favours the subtractive interaction channel forming lighter strings. 
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A comparison between the left and right panels of Fig. [21 shows that taking into account 
the kinematic constraints of equations ([H])-© is crucial for making quantitative predictions. 
By removing kinematically forbidden processes that would otherwise have been included, it 
introduces both an overall suppression factor with respect to self-interactions (which are not 
subject to these constraints) and a relative weighting among different cross-interactions. The 
overall suppression makes self-interactions more important, thus driving the average string 
densities closer to their self-interacting values (A^ ~ 14 for c = 0.23, ~ 980 for c = 0.023), 
while the relative weighting changes the abundances among different string species. The 
general trend that heavy, composite strings are systematically suppressed is observed, the 
dominant species being F-, D- and (±1, l)-strings. 

The corresponding scaling solutions for the network RMS velocities are shown in Fig. 
[31 Heavier, less populated strings are systematically slower, with the heaviest ones having 
negligible RMS velocities as their correlation length falls outside the horizon. This effect is 
exaggerated in these simplified models {b^j = 0), which ignore momentum transfers during 
binding processes; equations ([2]) for bfj = only take into account correlation-length-scale 
string curvature and Hubble damping, so, once a string species correlation length exits 
the horizon, Hubble damping dominates. At smaller scales, the velocity is controlled by 
short-scale structure, which is not captured in detail by these models. Allowing momentum 
transfers among different networks, introduces an extra source that can significantly affect 
the velocities of heavy strings, as will be demonstrated later. The basic feature that the 
heavier strings are systematically slower is consistent with the analytic analysis of the RMS 



velocities of junctions performed in 



34 



In Figs. [21 and [3l we have assumed that all self-interaction and cross-interaction processes 
have the same efficiency for all network components, that is, we have chosen a single value 
for all Cj's (the efficiency with which strings of type i self-interact producing loops) and 
a single value for all di/s (the efficiency with which strings of type i bind with strings 
of type j). This was only to demonstrate the general effects of cross-interactions and of 
junction constraints on the macroscopic properties of the network. In principle, the self- 
interaction and cross-interaction efficiencies depend on the types of strings involved and are 
related to the quantum intercommuting probabilities of each process, which can be estimated 



using string theory methods 



43 



44j |. These are strongly model-dependent, but the general 



picture resulting from such analysis is that intercommuting probabilities for F-strings are 
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FIG. 3: Velocity evolution for the string networks of Fig. [2](6^j = 0, gg = 0.3). F-strings are shown 
in black, D-strings in blue, (1, ±1) strings in red, (2, ±1) in green, (1, ±2) in orange, (3, ±1) in pink 
and (1, ±3) in yellow. Plots on the left do not take into account the effects of junction constraints, 



while those on right do; upper plots are for Cj = 0.23, dij = 1, and lower plots for Cj 



0.023, 



di 



0.1. 



likely to be in the range PpF £ [0.01, 0.1] while for D-strings Pdd ^ [0-1, 1]- For interactions 
between F- and {p, g)-strings, the corresponding probabilities are enhanced with respect to 
PpF by at least an inverse factor Qs, while for interactions between heavy composites one 
expects probabilities of order unity due to large multiplicities of the relevant Chan-Paton 
factors as well as due to a suppression of fluctuations for heavier strings. These different 
intercommuting probabilities must be reflected in the coefficients q, dij (although in a non- 
linear way ^^]) and can crucially affect the relative string abundances. In particular, the 
generic situation is Pff < -Pdd ^ cf < cd so these networks are generally dominated by the 
lightest (1,0) F-string components. It is worth therefore solving the system for somewhat 
more realistic values of the relevant parameters (Fig. Hj). The assumed values for the 
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parameters q and dij are shown in the legend. 
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FIG. 4: Evolution of normalised string density = pt'^ / fJ, in the radiation era, for a network 
with b^j = 0, but with a more realistic choice of the parameters Cj and dfj {gs = 0.3). Again, 

F-strings are shown in black, D-strings in blue, (l,ibl) strings in red, (2, ±1) in green, (l,ib2) 

1/3 

in orange, (3,ibl) in pink and (l,ib3) in yellow. Here, Cj = 0.23 x where, in the above 

order, Pj = (0.05,0.1,0.6,0.7,0.9,0.8, 1). The power 1/3 comes from the numerical simulations of 
Ref. [47]. The coefficients dij are: du = 0.4, di3 = 0.6, di^ = 0.65, dis = 0.8, dig = 0.75, diy = 
0.95, d23 = 0.5,^24 = 0.75, d25 = 0.85,^26 = 0.8,^27 = 0.95,^34 = 0.75,^35 = 0.87,^6 = 0.8,^37 = 
0.9,^45 = 0.9,^46 = 0.9,^47 = 0.95,(^56 = 0.9, d^j = 0.95 and dys = 0.95 (all unrelated are null). 
Plots on the right (left) do (do not) take into account the kinematic constraints on junctions. 



As already mentioned, the simple model described above is somewhat oversimplified in 
that, by switching off the interaction terms in the velocity equations ([2]), it neglects the 
effect of momentum transfer among different network components on the RMS velocity of 
each string type, which in turn affects string densities through the velocity dependencies 
of equations ([T]). Furthermore, the assumption that the energy differences associated with 
string binding can be efficiently radiated away, critically affects the cosmological evolution of 
the network. Indeed, if this assumption did not strictly hold, the simple system considered 
would artificially damp energy away leading to a possibly 'spurious' scaling regime. A more 
conservative approach would be to enforce energy conservation during string interactions by 
redistributing the relevant energy differences within the network through momentum transfer 
among interacting string segments. As mentioned earlier this is described by including the 
interaction terms in the velocity evolution equations ([2]) with b^j = d^y 
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In 



26| it was found that, for this system, scahng is not reached for generic choices of 



parameters. The energy hberated by binding processes remains within the network in the 
form of kinetic energy of newly formed segments, so in this case, and unhke 251], junction- 



forming processes alone can never lead to scaling behaviour. As in the case of single string 
networks, self-interactions and the related loop production that removes energy from the 
long-string network, are crucial for scaling to have a chance to develop. Nevertheless, scaling 
networks can be readily obtained in this case also, as long as the coefficients dij (binding 
processes) are somewhat smaller than the corresponding loop-production coefficients Cj. It 
is clear therefore that, in this case, including the effect of the constraints of section IIIII 
in the full system ([I])-([2]),(j4]),([5]), with b^j = d^j, can make the difference between network 
frustration and scaling. In particular, the constraints induce an extra suppression in the 
coefficients d'^j through the redefinition flTSl) . thereby providing a natural way to implement 
the requirement of d'^^ < q, as we saw in section [IVl a result which favours scaling. 

Fig. [5] shows the normalised string densities corresponding to Fig. [2] but this time for 
^ij = ^ij instead of b^j = 0. On the left panel (upper: q = 0.23, dij = 1; lower: Cj = 
0.023, dij = 0.1) string densities fail to reach scaling, as the energy liberated by binding 
processes, which exceeds the energy damped by loop-production, is redistributed within the 
network as kinetic energy. The suppression in the coefficients d^j = b^j due to the junction 
constraints readily solves this problem, leading to a scaling solution as can be seen in the 
right panels. Comparing the right panels of Figs. [5] and [2] it is clear that, after taking 
into account the kinematic constraints on string junctions, the effect of momentum transfer 
can be quantitatively significant, though the results are similar in both cases. The scaling 
network velocities corresponding to the right panel of Fig. [3] are shown in Fig. [61 Again, 
heavier, less populated species are slower, but as discussed above the effect is exaggerated 
in Fig. [31 

Finally, let us consider somewhat more realistic values for the coefficients q, dij, as in 
Fig. [H but now for = d^j instead of b^j = 0. This is shown in Fig. [71 The resulting 
scaling densities for the chosen intercommuting probabilities are very similar to those of 
Fig. [H The network is dominated by F-, D-, and (±1, l)-strings, in decreasing abundances, 
with a negligible contribution of heavier species. Note that the relative abundance of these 
three lightest species depends on the corresponding intercommuting probabilities, which are 
model-dependent. Overall, there is a significant enhancement in string density compared to 
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FIG. 5: Normalised string densities, as in Fig. [21 but now enforcing energy conservation at junction 
formation bf^ = d^^. On the left, where the effect of the constraints is not taken into account, the 
networks fail to scale as the junction formation terms dominate. Including the constraints (right) 
effectively suppresses these terms and scaling solutions are observed. 



field theory strings, whose normalised string density would scale at a value of ~ 14 in the 
radiation era. 

Now, all the plots in this section are for a string coupling Qs = 0.3. For smaller values 
of Qs the difference in tension between F- and D- strings becomes more pronounced and 
there is a larger energy gain associated with zipping. The coefficients favour even more 
strongly the 'unzipping' of heavier species into lighter ones, and the network becomes more 
strongly dominated by the lightest F-strings. For Qs = 1, on the other hand, F- and D- 
strings have the same tension and {p,q)-, (g,p)-strings are degenerate. There are then only 
four (instead of seven, in this truncation) non-degenerate tension species whose abundances 
again fall with increasing tension. Also, the plots produced here are for evolution in the 
radiation era. Network evolution in the matter era reaches scaling at lower string densities 
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FIG. 6: Network velocities corresponding to the scaling plots (right panel) of Fig. [5l 



(e.g. ~ 3 for c = 0.23, ~ 35 for c = 0.023, both in the self-interacting case d = 0). 



VI. DISCUSSION 



Understanding the evolution of a network of cosmic superstrings is a particularly impor- 
tant element in determining the cosmological properties of these intriguing objects. Through 
a combination of analytical and numerical approaches, evidence has been provided that such 
complicated networks which involve strings of many tensions which can form three-string 
junctions can reach scaling solutions, where the energy density in each type of string scales 
with the background energy density. In that regime, the lighter strings have a higher den- 
sity than the heavier strings. In this paper we have extended the semi analytic model of 



261] to include the kinetic constraints obtained in 33] for (p, g)-strings. These constraints 



have the effect of limiting the range of relative velocities and angles of approach which can 
lead to junction formation for interacting strings. Earlier analysis of the dynamics of string 
networks have failed to take these influential constraints into account, and in this paper we 
have established how important they are. In many ways, the nicest aspect of the analysis is 
that the net effect of the constraints can be accommodated by a rescaling of a single term 
in the one-scale evolution equations as seen in equation f[T5]) where the effective interaction 
strength between colliding strings is reduced. This has the knock on effect that the final 
scaling solution for the strings is increased beyond what would naively have been expected. 
Moreover, it also means that solutions that were previously found not to scale, now do en- 
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FIG. 7: Normalised string density (upper plots) and network velocities (lower plots) for the param- 
eters of Figs. [2] & [3] but for b^j = d^^, enforcing energy conservation at junction formation. Again, 
left (right) plots are without (with) the constraints. As was also the case in Fig. O these networks 
do not scale when the constraints are ignored (left). 



ter scaling regimes as seen in Figs. ^ and ([7]). The consequences of this for observations 
are important. An increased scaling density by a factor of 10 or so as seen in our analysis 
implies that there is more string than previously assumed. It is important that the con- 
sequences of this are followed up. We are currently determining the quantitative effects of 
these constraints on CMB anisotropies 49] . 
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